A model displaying extremely inhomogeneous matter 
distribution in General Relativity 
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We consider a toy metric in four dimensional space-time defined in 
terms of a recursive hierarcliical prescription. The matter distribution 
turns out to be extremely inhomogeneous. Surprisingly, for very large 
tJ- I samples the average mass density tends (very slowly) to a constant. 

' There is no trace of fractal dimension left. 
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! Although most authors tend to believe that the large scale distribution 

^ ' of matter should be constant as predicted by Einstein's cosmological prin- 

ciple [Tj there is also a belief that such an idea should be further checked 
p • against observation j2]. Some authors argued quite convincingly that the 

D ! large scale distribution of matter should be fractal, quite the oposite to that 

^ ' expected by the cosmological principle. Inhomogeneous cosmological mod- 

els have been searched by many authors |S] 111 [HI [Z| [H]- The fact that 
^ ■ detailed maps show quite strong inhomogeneities has been a motivation for 

. long controversies P [TUj . 

Here, we construct a metric defined in terms of a hierarchy, building 
an expectation that the matter distribution obtained from the corresponding 
Einstein's Equations are extremely inhomogeneous, contradicting the cosmo- 
logical principle, and probably presenting fractality properties, with a non 
trivial effective dimension. As it turns out, some of these conclusions are 
false, as we shall see. 

Our starting point is the diagonal metric defined in terms of a distribution 
of points. A point in space-time is only defined at the integers. Indeed, we 
expect that a metric describing a fractal distribution cannot be continuous 
at any point, therefore it can only be defined on a discrete set. With this 
motivation in mind, we define, in a given direction i (either 1, 2 or 3) the 
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distance between two points as pictorically represented in Fig. 1, with the 
metric 
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where the metric is defined on all integers, depending on their decomposition 
in terms of powers of 2, as below, 

g-^^[x) = a^^ , with X = 2^=+^ + 2^= - 1 

922{y) = with y = 2'+i + 2'-l (2) 

^33 (x) = a^"* , with z = 2"^+^ + 2™ - 1 . 

Above, I and m are nonnegative integers and the metric is defined as unit 
when the entry is an even number. This defines the metric (see Fig. 1). 

Computing the Christoffel symbols and subsequently the curvature tensor 
requires some care, since we are not dealing with derivatives of functions, 
but differences of functions defined on a discrete space. The result for the 
Christoffel symbols is 

r^o = -§, r^2 = ^^(i-i?") , (3) 

a at 1 
_mda 3 _ (-1)% ^ _ 33. 
adt' ^33- 2 ^ > ■ 

A naive construction of the Riemann tensor may lead to a result incompatible 
with the general symmetry properties required for that tensor ^1]. We thus 
define it as 

R = 1 r d'^9\u d'^9iJ.iy d'^9KX d'^9tiK -, r-pr? -po- _-p»? -po- t 

which in the continuous version corresponds exactly to the usual definition, 
and here has the correct symmetry properties with respect to the interchange 
of indices. The Ricci tensor is readily obtained, with the result (the non 
diagonal components vanish) 

Rm = ^(^ + ^ + "^)^ + ^(^) [kik-l) + l{l-l) + mim-l)] , 
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Figure 1: Diagram showing the definition of a hierarchical metric as a func- 
tion of the distance. In Eq. |21we defined a metric which reahzes the above 
properties for nearby points. The arc segment with a power of a shown here 
represents the distance between any element from one side of the arc and 
any element from the other side. 
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[kl + l{l - I) + Im]} 
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m 1 ( da\ . . , , . , 
R33 = ^33|— + ^| — I [mk + ml + m{m - 



The matter density can be now readily calculated by means of the Einstein 
Equations, and we obtain 

1 / 1 da\ ^ 

Too = p=——\-—] (kl + lm + mk)=po(kl + lm + mk) . (5) 
SvrG \a at J 

The result shows severe inhomogeneities. A graph showing p{x, y, 0) exem- 
plifies the situation in Fig. 2. 

We further questioned about the fact of whether such an inhomogeneous 
distribution could be described by a fractal. We numerically computed the 
limit 

l.m^ , (6) 

r— »oo fO- ^ 

where N{r) = J2o<x,y,z<r Pi^iUi ^) finding the value of d where the limit 
exists. This is the usual definition of a fractal dimension 12j. Some results 
are sumarized in Figure 3. 

It turned out, with some surprise, that the effective dimension of the 
distribution is (i = 3, that is, the hierarchical metric was not enough to lead 
to a hierarchical distribution of matter. Nonetheless, the matter distribution 
is very inhomogeneous, and does not seem to homogenize very soon. Whether 
such a simplified model can have any relation to the real matter distribution 
of the universe is not questioned, but the stability of the dimensionality 
of the effective distribution of matter is quite remarkable. This may give 
support to models displaying highly inhomogeneous behaviour either in the 
nonrelativistic domain, such as [6^, or up to a certain cosmological distance, 
but smooth out their behaviour at large enough scales, as predicted by the 
cosmological principle, very slowly approaching the canonical value d = 3 for 
the effective dimension. 
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Figure 2: Matter distribution for a fixed value of the ^-component, displaying 
the observed inhomogeneity. 
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Figure 3: Evolution of the average density as further data are included, 
according to the limit procedure reflimit. The lower diagram corresponds to 
d = 3.25, the next one to d = 3.05, and the highest one to d = 3.00. We note 
that the highest one corresponds to a constant, thus giving the correct value 
of d. We also note the fact that we need a long time to achieve a cpnstant 
value for the limit. 
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